Draft version February 14, 2008 

Preprint typeset using I£T^X style cmulateapj v. 10/09/06 



oo 
O 
O 
(N 

X> 

IX, 
(N 



43 

6 



> 
m 

(N 
O 
OO 
O 



X 



SPITZER IRS OBSERVATIONS OF DISKS AROUND BROWN DWARFS IN THE TW HYDRA ASSOCIATION 1 

A. L. Morrow 2 , K. L. Luhman 2,3 , C. Espaillat 4 , P. D'Alessio 5 , L. Adame 6 , N. Calvet 4 , W. J. Forrest 7 , B. 
Sargent 7 , L. Hartmann 4 , D. M. Watson 7 , & C. J. Bohac 7 

Draft version February 14, 2008 

ABSTRACT 

Using SpeX at the NASA Infrared Telescope Facility and the Spitzer Infrared Spectrograph, we have 
obtained infrared spectra from 0.7 to 30 /im for three young brown dwarfs in the TW Hydra Association 
(r ~ 10 Myr), 2MASSW J1207334-393254, 2MASSW J1139511-315921, and SSSPM J1102-3431. The 
spectral energy distribution for 2MASSW J1139511-315921 is consistent with a stellar photosphere for 
the entire wavelength range of our data while the other two objects exhibit significant excess emission 
at A > 5 fjjn. We are able to reproduce the excess emission from each brown dwarf using our models 
of irradiated accretion disks. According to our model fits, both disks have experienced a high degree 
of dust settling. We also find that silicate emission at 10 and 20 fiva is absent from the spectra of 
these disks, indicating that grains in the upper disk layers have grown to sizes larger than ~ 5 ^m. 
Both of these characteristics are consistent with previous observations of decreasing silicate emission 
with lower stellar masses and older ages. These trends suggest that either 1) the growth of dust 
grains, and perhaps planetesimal formation, occurs faster in disks around brown dwarfs than in disks 
around stars, or 2) the radii of the mid-IR-emitting regions of disks are smaller for brown dwarfs than 
for stars, and grains grow faster at smaller disk radii. Finally, we note the possible detection of an 
unexplained emission feature near 14 /im in the spectra of both of the disk-bearing brown dwarfs. 
Subject headings: accretion disks - planetary systems: protoplanetary disks - stars: formation - 
stars: low-mass, brown dwarfs — stars: pre-main sequence 



1. INTRODUCTION 

One of the first steps toward the formation of planets is 
the growth and settling of dust grains in a cir cumstellar 
accretion disk ( Wciden schilling fc Cuzzi|["l993h . A com- 
mon diagnostic of gra in evolution in dis ks is the silicate 
emission near 10 /im (jNatta et a,l.ll2000h . Measurements 
of this feature have been performed for young stars (r ~ 
1 Myr) acr oss a wide range of spectral types, from early B 
to late M dFurlan et al1l2005L iKes sler-Silacci et~aD l2006t 
ISicilia-Aguilar et al.ll2007h ~In these data, silicate emis- 
sion becomes weaker at later spectral types, which could 
be explained by more advanced grain evolution or sedi- 
mentation in disks around low-mass objects or a lumi- 
nosity dependence of the ra dius in the disk at which 
silicate emission is produc ed (|Kessler-Silacci et al.ll2007t 
ISicilia-Aguilar et al.ll200l . 

Recent studies also have begun to explore the evo- 
lution of silicate emission with time for low-mass sys- 
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terns. Observations with the Spitzer Space Telescope 
have found that disks around low-mass stars and brown 
dwarfs exhibi t weaker silicate em ission in Upper Scorpius 
(r ~ 5 Mvr. iScholz et aLll2007f) than in Chamaeleon I 



r ~ 3 Myr lApai et al.ll2005h and Taurus (r ~ 1 Myr, 



iFurlan et a l. 2005), which is consistent with the growth 
of grains to larger sizes, the settling of dust to midplanc, 
or both, as time goes on. However, according to a recent 
study, this correlation between silicate emission and age 
does not apply to the brown dwa rf 2MASSW J1207334- 
393254 (henceforth 2M 1 207-3932 .lGl^ l2002h in th e TW 
Hya Association (TWA, r ~ 10 Mvr. iMamaiekl 120051 : 
iBarrado v Navascues etall l2006t ) . iRiaz fc GiziT(|2007l ) 
constructed a mid-infrared (IR) spectral energy distribu- 
tion (SEP) for this objec t usin g broad-band pho tometry 
from ISterzik et all (|2004D and lRiaz et"afl (|2006l) . which 
seemed to indicate the presence of silicate emission. As 
a result, they concluded that the disk around 2M 1207- 
3932 has experienced little dust processing relative to 
stars at the same age. 

To investigate the evolution of silicate emission in 
brown dwarf disks more definitively, we have performed 
mid-IR spectroscopy on 2M 1207-3932 and two other 
brown dwarfs i n TWA with the S pitzer Infrared Spec- 
trograph (IRS: iHouck et alj [2004T ). In this Letter, we 
describe these observations and supporting near-IR spec- 
troscopy (§ [2]), examine these data for silicate and con- 
tinuum emission from circumstellar dust (§ I3.ip . and fit 
these data with the predictions of disk models 13. 2|) . 

2. OBSERVATIONS 

2.1. Near-infrared Spectroscopy 

We obtained low-resolution near-IR spectra of 
2M 1207-3932 and two other brown dwar fs in TWA 
2MASSW J1139511-315921 (2M 1139-3159, [Giiii[2002) 
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and S SSPM J1102- 3431 (SS 1102-343 1, iScholz et all 
l2005aj) . with SpeX (jRavner et all l2003f ) at the NASA 
Infrared Telescope Facility (IRTF) on the night of 2005 
December 14. These data were re duced with the Spex- 
tool package (Cushing ct al. 2004) and corrected for tel- 
luric absorption (|Vacca et al.l 120031 ). The final spectra 
extend from 0.8-2.5 /jm and exhibit a resolving power of 
R = 100. We flux calibrated the spectra using photome- 
try at J, H, and K s from the Point Sourc e Catalog of the 
Two- Micron All-Sky Survey (2MASS, ISkrutskie et~"aTl 
2006). To measure spectral types from these data, we 
compared them to SpeX data for young objects that have 
been classified at optical wavelengths. The strengths 
of the TiO, VO, and H2O absorption bands indicate 
a spectral type of M8.5 for each object, which is con- 
sistent with the optical types of M8- M8.5 reported by 
iGizisI (|2002[ ) and IScholz et all (|2005al) . The spectra of 
2M 1139-3159 and SS 1102-3431 are slightly redder than 
the spectrum of 2M 1207-3932 in a manner that is consis- 
tent with reddening by interstellar dust. Therefore, we 
dereddened the spectra of 2M 1139-3159 and SS 1102- 
3431 by Ay — 0.8 and Ay = 0.5, respectively, to match 
the data for 2M 1207-3932. 

2.2. Mid-infrared Spectroscopy 

We obtained low-resolution mid-IR spectra of 2M 1139- 
3159, 2M 1207-3932, and SS 1102-343 on 3 July 2005, 
29 July 2006, and 3 July 2005, respectively, with IRS 
aboard Spitzer as a part of the Guaranteed Time Obser- 
vations of the IRS instrument team. We used both low- 
resolution IRS modules, Short-Low and Long-Low, which 
cover 5.3-14 and 14-40 /im, respectively, with a resolution 
of A/AA ~ 90. The total exposure time for each target 
was ~ 4000 sec. The spectra were extracted and cali- 
brated from the basic calibrated data usin g the standard 
SMA RT data analysis package for IRS (Higdo rTet al.l 
120041 ). The spectra were reduced with the methods that 
have been previousl y applied to IRS d ata for low-mass 
members of Taurus (|Furlan et al.ll2005T ). 

3. ANALYSIS 

3.1. Disk Emission 

In Figure □ we plot the SED of 2M 1139-3159 from 
0.7 to 24 /xm usin g photometry from 2MASS and Spitzer 
(R iaz et al.ll2006| ) and the spectra that we obtained with 
SpeX and IRS. For comparison, we include the SED of 
a field dwarf with a similar spectral type (VB 10, M8V) 
after scaling it to match 2M 1139-3159 at J, H, and 
K s . 2M 1139-3159 is slightly brighter that the field 
dwarf beyond 3 /im, but this is probably a reflection of a 
small difference in the photospheric near- to mid-IR col- 
ors between pre-main-sequence objects and field dwarfs 
rather than excess emission from dust given that the 
mid-IR slopes are the same for the two objects. There- 
fore, we believe that the SED of 2M 1139-3159 repre- 
sents a good estimate of the SEDs of the stellar photo- 
spheres of the other two brown dwarfs, 2M 1207-3932 and 
SS 1102-3431. We compare the SpeX and IRS spectra 
of these three objects in Figure [2 Relative to the pho- 
tospheric SED of 2M 1139-3159, the SEDs of 2M 1207- 
3932 and SS 1102-3431 exhibit significant excess emis- 
sion longward of 5 /xm. These results for 2M 1207-3932 
and 2M 1139-3159 are consistent with those of previ- 
ous studies, which have found evidence of a disk for 



2M 1207-393 2 but not for 2M 1139-3159 based on Spitzer 
photometry dRiaz et al.ll2006D. ground-based ph otometry 
(Pavawardhana et al.ll2003t ISterzik et al.ll2004l ). and Ha 
emission (jMohantv et al.ll2003ll2005f) . For SS 1102-3431, 
the only previous constr aint on the p r esence of a disk is 
the Ha spectroscopy bv IScholz et al.l <|2005bt ). who con- 
cluded that accretion is occurring at a very low level, 
if at all. Thus, our IRS observations provide the first 
definitive detection of a disk around this object. 

In addition to detecting the presence of excess emis- 
sion, previous mid-IR photometric measurements of 
2M 1207-3932 have been used to constrain the strength of 
the silicate emission feature at 10 [im. Bas ed on ground- 
based photometry at 8.7 and 10.4 /im, ISterzik et al.l 
(|2004f) found that 2M 1207-3932 did not exhibit strong 
silicate emission. However, 2M 1207-3932 was slightly 
brighter at 10.4 /im than at 8.7 /im, which was in- 
terpreted as evide nce of modest silicate emission by 
iRiaz fe Gizisl ((2007). As shown in Figure[2j silicate emis- 
sion at 10 and 20 /im is absent in the IRS spectrum 
of 2M 1207-3932. The same is true for the other disk- 
bearing brown dwarf in our sample, SS 1102-3431. In 
contrast to these brown dwarf disks, the disks around 
stars in TWA do produc e silicate emission (|Uchida et al.l 
12003 Ionian et al.ll2007h . 

Although silicate emission at 10 and 20 /im is not 
present in the IRS data for 2M 1207-3932 and SS 1102- 
3431, it appears that a weak emission feature may be 
detected near 14 /im in the spectra of both objects. Be- 
cause of its large width, the feature is probably produced 
by a mineral rather than a gas. We cannot identify a 
plausible source of this emission that is consistent with 
the available data for these disks (e.g., large grain sizes, 
absence of other emission lines). 

3.2. Disk Model 

We have modeled the mid-IR excess emission from 
2M 1207-3932 and SS 1102-3431 as arising from irradi- 
ated accretion disks by fol l owing the procedures from 
iD'Alessio et "all fflQ98L Il999l l200lT I2006T) . For the stellar 
photosphere of each object, we have adopted an effective 
temperature of 2555 K, which is estimated by combin- 
ing a spect ral type of M8. 5 (jj I2.ip with the temperature 
scale from iLuhman et alJ ((2003) . A brown dwarf with 
this temperature at the age of the TW Hya association 
(10 Myr) is predicted to have a bolometric luminosity of 
0.0024 Lq and a mass of .025 M q by the evolutionary 
models bv IChabrier et ail ((2000). These values of tem- 
perature and luminosity correspond to a stellar radius of 
0.25 i?©. For the disk calculations, we adopted a uni- 
form a ccretion rate of M = 10~ n M Q yr _1 (Scho lz et alJ 
2005b) . W e include dust settling i n the disk following the 
methods of ID'Alessio et all (|2006l ). In short, two popula- 
tions of grains co-exist in the disk, both with size distri- 
butions given by n(a) ~ q~ 3,5 where a is the grain radius 
(|Mathis. Rumpl. fc Nordsiecklll977D . The minimum ra- 
dius for both populations is 0.005 /im; grains around the 
midplane have a maximum radius of a max = 1 mm, while 
a m ax in the upper layers is adjusted to produce the best 
fit to the SED. In addition, the dust-to-gas mass ratio 
of the small grains is parameterized in terms of e, which 
is the ratio normalized by the standard interstellar value 
of ~ 0.01. The model includes an inner disk wall illu- 
minated by the central brown dwarf and located at the 



Brown Dwarf Disks in TW Hya 



3 



dust sublimation radius of 3.3 R*. The outer radii of the 
disks around 2M 1207-3932 and SS 1102-3431 are not 
constrained by the available data; we adopt a value of 
Rout = 50 AU for each disk. 

We have calculated SEDs for a range of values of the 
viscosity parameter (a), e, disk inclination (i), and wall 
height (z wa u) and compared the results to the observed 
SEDs of 2M 1207-3932 and SS 1102-3431. The SEDs 
of the best model fits are shown in Figs. [3] and 2] The 
data for both objects are well-matched by a — 0.001, 
e = 0.001, and i — 60°. We are able to reproduce the 
small difference in the SEDs of the two objects (Figure[2]) 
by using z waU = 1 if for 2M 1207-3932 and z waU = 2 H 
for SS 1102-3431, where H is the disk scale height, which 
is 1.15 x 10~ 4 AU for the adopted sublimation temper- 
ature of 1400 K. These differences in z wa u may indicate 
that the inner disks have different degrees of dust set- 
tling. We note that other combinations of values for 
a, e, and i also produce reasonable fits to the observed 
SEDs. Additional data, such as photometry at longer 
wavelengths, are needed to further constrain these pa- 
rameters. However, the relatively blue slopes of the mid- 
IR SEDs of 2M 1207-3932 and SS 1102-3431 definitely 
indicate a large degree of dust sett ling to the disk mid - 
plane (i.e., small e in our treatment. IWatson et al.ll2008f ). 

The absence of silicate emission (Figure [2]) indicates 
that grai ns have grown signific antly in the upper disk 
layers (cf. iD'Alessio et al.|[2006t) . We show in Figured a 
series of disk models in which we have varied a max in the 
upper disk layers from 0.25 /im (ISM grains) to 100 /im. 
As expected, the silicate emission becomes weaker as 
a ma x grows. Based on the comparison to our models, the 
absence of silicate emission in the IRS data for 2M 1207- 
3932 and SS 1102-3431 indicates that grains in the upper 
disk layers have grown to sizes of a max > 5 /im. 
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4. CONCLUSIONS 

We have presented near- and mid-IR spectroscopy for 
three young brown dwarfs in the TW Hya association. 
Two of these objects, 2M 1207-3932 and SS 1102-3431, 
exhibit significant mid-IR emission above that expected 
from a stellar photosphere. We have successfully repro- 
duced the excess emission from each brown dwarf with 
an irradiated accretion disk model that includes dust set- 
tling. Both disks exhibit high degrees of dust settling to 
the midplane based on the relatively blue mid-IR slopes 
of their SEDs. In addition, our IRS spectra reveal an ab- 
sence of silicate emission at 10 and 20 /im in both objects, 
which indicates that the disks have experienced signifi- 
cant grain growth in their upper layers (a max > 5 /im). 
These results for 10-Myr-old brown dwarfs support and 
extend the previously observed trend of decreasing sili- 
cate emission with lower stellar masses and older ages. 
This trend may indicate that grain growth and planetesi- 
mal formation occur more rapidly in disks around brown 
dwarfs than in disks around stars, or that grains grow 
faster at smaller disk radii, and it is at smaller radii where 
mid-IR emission is produced for objects at lower mas ses 
(|Kessler-Silacci et al J 120071: ISicilia-Aguilar et al.ll2007t ). 
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Fig. 1. — SEDs of the young brown dwarf 2M 1139-3159 (lines and filled circles) and the field dwarf VB 10 (open circles). The SED 
of 2M 1139-3159 is consistent with that of the field dwarf and does not show excess emission at long wavelengths from a circumstellar 
disk. The spectra of 2M 1139-3159 were obtained with SpeX and IRS in this work. The photometric measurements are from 2MASS, 
IPatten et all ll200rT) , and lRiaz etTaTI ([2003) . The SED for VB 10 has been scaled to match that of 2M 1139-3159 at J, H, and K s . 
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Fig. 2.— SEDs of the young brown dwarfs SS 1102-3431, 2M 1207-3932, and 2M 1139-3159 (red, black, and blue lines). Relative to 
2M 1139-3159, SS 1102-3431 and 2M 1207-3932 exhibit excess emission at long wavelengths, indicating the presence of circumstellar disks. 
The difference between the SEDs of SS 1102-3431 and 2M 1207-3932 can be explained as a difference in the heights of their inner disk walls 
(Figs. |3] and gj. The SEDs of SS 1102-3431 and 2M 1139-3159 have been scaled to match that of 2M 1207-3932 at 1-2.5 (an. 
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Fig. 4. — SED of 2M 1207-3932 compared to a model for its circumstellar disk. The model parameters are the same as for SS 1102-3431 
in Figure \3\ except that z wa u = 1 H and four values of a ma x are shown (0.25, 1, 5, and 100 fim). The predicted strength of the silicate 
emission near 10 /an is weaker for larger values of ama x- We adopted the SED of 2M 1139-3 159 to represent the stellar photosphere of 
2M 1207-3932. The photometric measurements are from lSterzik et al.| iffijoU, open circles) and[Riaz et al. (2006], filled circles). 



